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Abstract 



We introduce the notion of characteristic function of a quaternionic 
matrix, whose roots are the left eigenvalues. We prove that for all 
2x2 matrices and for 3x3 matrices having some zero entry outside 
the diagonal there is a characteristic function which is a polynomial. 
For the other 3x3 matrices the characteristic function is a rational 
fT^ . function with one point of discontinuity. We prove that Hamilton- 

^ ' Cayley theorem holds in all cases. 

t^ ' Keywords: quaternion, left eigenvalue, characteristic function, Hamil- 

ton - Cayley theorem 
MSC: 15A33, 15A18 



1 Introduction 



Very little is known about left eigenvalues of n x n quaternionic matrices. 
rS I F. Zhang's papers |9l|10j review their main properties as well as some patho- 

C^ • logical examples, see also [3]. For n = 2 the explicit computation of the left 

spectrum is due to L. Huang and W. So [4], while the authors studied the 
symplectic group in [SJ |6] . 

In 1985, R. M. W. Wood [8] proved, by using homotopic methods, that 
every quaternionic matrix has at least one left eigenvalue. At the end of his 



paper. Wood notes that "in the 2x2 case of the matrix 



a 



b 



there is a 
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partially defined determinant b — ac ^d and partially defined characteristic 
equation 

Ac-^A - Xc"^d - ac'^X -b+ [ac~'^d] = (1) 

which reduces the eigenvalue problem to the fundamental theorem [of alge- 
bra]. The difficulties start with 3x3 matrices". 

In this paper we introduce a definition of characteristic function for a 
quaternionic matrix, which generalizes the usual characteristic polynomial in 
the real and complex setting. In particular, its roots are the left eigenvalues. 
Explicitly, we say that /li : BI — > IHI is a characteristic function of the matrix 
A & Ai{n, H) if, up to a constant, its norm verifies that |/i(A) | = Sdet(^— A/) 
for all A € H, where Sdet: A4{n,M) — )• [0, +oo) is Study's determinant. As 
we shall see, this definition fits naturally with Equation ([1]), as well as with 
the method proposed by W. So in [^ to compute the left eigenvalues when 
n = 3. 

Then we discuss Hamilton-Cayley theorem in this setting. Our main 
result is as follows. 

Theorem A. For any quaternionic matrix A E A4(n,M), n < 3, there 
exists a characteristic function fj, whose extension to a map /i: A4{n,M) — )• 
Ai{n,M) verifies Hamilton-Cayley, that is /u(^) = 0. 

For n = 2, a characteristic function like that in ([T|) is a polynomial /u(A) 
for which it is easy to check that /u(^) = 0. It follows that 

Ac'^A = Ac~^d + ac'^A + (6 - ac~'^d)I, 

which generalizes the well known formula A^ = {ticA)A — (det A)I in the 
commutative setting. When n = 3 and the matrix has some zero entry 
outside the diagonal, we shall find a polynomial characteristic function that 
verifies Hamilton-Cayley. Otherwise, there is a characteristic function which 
is, outside a point of discontinuity, a rational function. We are able to 
extend it to a map fi: M.{n,M.) — )• M.{n,M) and we prove by brute force 
that Hamilton-Cayley is verified too. 

At the end of the paper we discuss another possible definition of charac- 
teristic function. 

2 Preliminaries 

We consider the quaternionic space H" as a right vector space over H. Two 
square matrices A,B£ Ai{n,M.) are similar if B = PAP~^ for some invert- 
ible square matrix P. 



If ^ is a quaternionic n x n matrix, let us write A = X + jY , with 
X,Y £M{n,C), and let 



ciA) 



X -Y 
Y X 



eX(2n,C) 



be its complex form. We have c{A ■ B) = c{A) ■ c{B) , c{A + B) = c{A) + c{B) 
and c{tA) = tc{A) if t G M. In particular, A is invertible if and only if c{A) 
is invertible. Moreover, det c{A) > is a nonnegative real number, so we 
can define the Study 's determ,inant of A as 

Sdet(A) = (detc(^))^/2 >o. (2) 

For complex matrices, Sdet equals the norm of the complex determinant, see 
[U 12] for a general discussion of quaternionic determinants. The following 
properties are immediate: 

1. Sdet(^ • B) = Sdet{A) ■ Sdet(B); 

2. A is invertible if and only if Sdet(j4) > 0; 

3. if A,B are similar matrices then Sdet(A) = Sdet(-B). 

We also need the following result. 

Lemma 2.1. For a matrix with boxes M,N of size m x m and n x n 
respectively we have 



Sdet 



M' 

N * 



Sdet(M) • Sdet(iV). 



It follows that Sdet(A) = \qi • • • Qnl when yl is a triangular matrix, with 
qi, . . . ,qn being the elements of the diagonal. 

Sometimes we shall permute two columns or rows of the matrix A. Or 
we shall add to a column a right linear combination of the columns. This 
will not affect the value of the determinant because the matrices of the type 



P 



"0 1 0" 




"1 0' 


1 


or P = 


a 1 


1 




13 1 



verify Sdet(P) = 1. 



Remark 2.2. Up to the exponent 1/2 in ([2]), this is the same determinant 
that the one in Theorem 8.1 of [9] that we shall refer to later in Sect. [71 The 
exponent is normalized in order to have Sdet(^) = \qi- ■ ■ qn\ for a diagonal 
matrix A = diag((7i, . . . , qn). 



3 Left eigenvalues and characteristic functions 

A quaternion A G H is said to be a left eigenvalue of the matrix A £ A4{n, M) 
if Av = Xv for some vector v G H", v ^ 0. Equivalently, the matrix A — \I is 
not invertible, that is Sdet(^ — XI) = 0, where Sdet is Study's determinant 
defined in Section [2j 

Definition 3.1. A map /i : H — > H is a characteristic function of the matrix 
A G M{n,M) if, up to a constant, |/u(A)| = Sdet(^ - XI) for all A G H. 

Notice that A is a left eigenvalue of A if and only if /i(A) = 0. 

Remark 3.2. It is well known that the left spectrum is not invariant under 
similarity. However, if P is a real invertible matrix then Sdet{PAP^^—XI) = 
Sdet(^ — XI), so A and PAP^^ have the same characteristic functions. 

Example 3.3. Diagonal and triangular matrices. 

If A = diag((7i , . . . ,qn) then /u( A) = ((jr„ — A) • • • (gi — A) is a characteristic 
function. Analogously for triangular matrices. 

Example 3.4. 2x2 matrices. 



Let A 



a b 
c d 



G M{2,m). If 6 = then Sdet(^) = \da\ and the map 



^(A) = {d — X){a — X) is a characteristic function. If 6 ^ we have 

A 



b 

c — db~^a d 



so 

Sdet(^) = \b\\c-db-^a\. 

Consequently, we consider the characteristic function 

^(A) = c-(d-A)6-^(a-A). (3) 

Obviously, the characteristic function of a matrix is not unique. For 
instance, by permuting rows and columns we can obtain /i(A) = b — (a — 
X)c~^{d — X) which is Wood's function in Equation ([T]) (there is a misprint 
in the original article). However, as we shall see in Section HI it is preferable 
to take minors starting from the top right corner as we do. 



4 Characteristic function of 3 x 3 matrices 

a b c 
Now let ^ = / g h be a 3 x 3 quaternionic matrix. The computation of 

p q r 
Sdet(^) can be done as follows (a similar algorithm is valid for any n > 3). 

4.1 Caser2 = 3,c^0 

First we consider the generic case when c ^ 0. In this case we can create 
zeroes in the first row, 



A 



c 

f — hc^^a g — hc^^b h 
p — rc~^a q — rc~^b r 



By Lemma |2. II and the 2x2 case, it follows: 
Proposition 4.1. If c y^ 0, then Sdet(yl) is given: 

1. when g — hc^^b / 0, by 

\c\ ■ \g — he" b\ ■\p — rc~ a — {q — rc~ b){g — hc~ b)~ (/ — hc~ a)|; 

2. when g — hc^^b = 0, by 

\c\ ■ \q — rc~ b\ ■ \f — hc~ a\. 



Corollary 4.2. Let us call \q = g — he ^b the pole of A. Then 

Sdet{A - XqI) = \c\ -Iq-ir- Xo)c^^b\ ■ \f - hc'^ia - Ao)|. (4) 

By applying Prop. 14.11 and Cor. 14.21 to A — XI we find the following 
characteristic function of A. 

Definition 4.3. When c 7^ 0, a characteristic function for the 3x3 matrix 
A can be defined as follows: 

1. a Xq = g — hc~^b is the pole of A, 

/^(Ao) = {q-{r- Xo)c-'b) (/ - hc-'ia - Aq)) ; 



2. otherwise, 

fi{X) = (Ao-A)[(p-(r-A)c-i(a-A))- 

{q-(r- X)c~'b) (Ao - A)-i (/ - hc-\a - A))] . 

Remark 4.4. In [7], W. So proved that the left eigenvalues of a 3 x 3 matrix 
can be computed as roots of certain polynomials of degree < 3. Even though 
our computation is different from his, we obtain that the function in Def. 14.31 
is exactly So's formula in [7, p. 563]. This is why we have chosen to compute 
determinants starting from the top right corner. 

4.2 Caseri = 3,c = 

We briefly review what happens when c = 0. First, if both b,h = we have 
a triangular matrix, then we can take 

MA) = (r-A)(ff-A)(a-A). (5) 

If 6 = but h ^ we can reduce to the 2x2 case by Lemma [2. H so we take 

MA) = (<7-(r-A)/i-i(<7-A))(a-A). (6) 

Finally, if 6 7^ we can (see the proof of Theorem 16. 3p create a zero in the 
left top corner of ^ — AI and then permute the second and last column, in 
order to reduce the matrix (^4 — XI)P to the 2x2 case. Alternatively, we 
can simply permute the second and last column and the second and last 
row of ^, in order to obtain a matrix PAP~^ with the same characteristic 
function, to which Subsection 14.11 applies. Notice however that with the 
latter method we obtain a rational function, not a polynomial. 

5 Continuity 

The following example shows that the characteristic function fj, in Definition 
14.31 mav not be continuous, even if its norm |//| is a continuous map. 

Let 

"0 i 1 

A= 3i-k 1 

_ k -1 + j + k 0_ 

Its pole (see Cor. 14. 2p is Aq = —i and 

^(Ao) = (j + k){2i -k) = l-i + 2j -2k. 



However, for A 7^ Aq we have 

^(A) = {-i - A) (A: - A^ - (-1 + j + A; + Xi){-i - X)-\3i - k + X)) , 
and by taking X = —i + ej, e G M, with e — ?> 0, we obtain 
hm ii{-i + ej) = l + i + 2j + 2k ^ fi{Xo). 

In fact, the hmit 

hm fi{—i + eq) = —q(j + k)q'' (2i — k) 

depends on q, so hniA^Ao f^W does not exist. 

It is an open question whether it is always possible to find a continuous 
characteristic function. 

6 Hamilton- Cayley theorem 

We now discuss Hamilton-Cayley theorem. 



6.1 Case n = 2 
Theorem 6.1. Let A 



a b 
c d 



be a 2 X 2 quaternionic matrix. Let /u(A) = 

c — {d — X)b^^{a — A) be the characteristic function defined in ^. Then 
fi{A) = 0. 



Proof. We have 

c 0' 
c 



d — a —b 
-c 



b'^ ■ 
b-^ 



-b ■ 




"0 0" 


-c a — d 








D 



Corollary 6.2. Ab-'^A = Ab'^a + db'^A + (c - db-^a)L. 

6.2 Case n = 3, c = 

For n = 3, a direct computation will show that Hamilton-Cayley theorem is 
true when c = (see Section H]) . 



Proposition 6.3. Let A 



a b 
f 9 h 
p q r 
tion defined in Subsection\4.^ Then fJ.{A) = 0. 



. Let fi{X) be the characteristic func- 



Proof. If 6, /i = we take formula ([5]), so fJ-{A) equals 



r — a 







g-a 
























-f 


r-g 




-/ 




-f a-g 


= 








. -P 


-q 0_ 




_ -P -q 9-r_ 




— p — g a — r 










If 6 = 


= 0, /i 7^ we take formula i^, then we check 






r — a 


" 




'g-a ' 


" 









-f 


r-g -h 


h-' 


-f -h 


-/ a-£ 


1 -h 




= 


. -P 


-q 0_ 




. -P -q i9-r)_ 


_-p -g 


a — r 










" 











q 


-/ a-g 


-h 


1 


















-p -q 


a — 


r 







that is, (r/ - A)h~^{gl - A){al - A) = q{al - A), hence n{A) = 0. 
If 6 7^ 0, we have 



6" 

f - {g - \)h-\a - X) h g-X 
p — qb~^{a — X) r — X q 



Sdet{A - XI) = Sdet 



so we are in the 2x2 situation (see Lemma l2. ID . First, assume h = and 
let us take ^(A) = {r - X) {f - {g - X)b~^{a - A)). We check 



r — a —b 
-/ r-g 
-p -q 



g — a —b 
-/ 

-p -q g-r 



■ -b ■ 

-f a-g 
—p —q a — r 

r — a —b 
-/ r-5 /, 
—p —q 

that is {rl - A)(gl - 71)6"^ (a/ - A) = {rl - A)f, hence fj.{A) = 0. 
On the other hand, if /i 7^ we take 

m(A) =p- qb-\a - A) - (r - X)h~^ [f - {g - X)b-\a - A)) . (7) 

Then we compute 

pI-qb-^{aI-A)-{rI-A)h-^f = 

^p-{r-a)h-^f q + bh-^f 

qb^f + fh-^f p-qb-\a-g)-{r-g)h-^f qb-^h + f 



qb ^p + ph ^f 



qb ^q + qh ^ f 



P — qb ^{a — r) 



and we check it equals 

r — a —b 
-/ r - g -h 
—p —q 

hence /i(A) = 0. 



h-' 



g — a —b 
-/ -h 
-P -q 9 -r 



0-6 " 

-f a-g -h 
—p —q a — r 



-{rl - A)h'\gl - A)b~\al - A), 



D 



Lemma 6.4. Let A be a quaternionic matrix such that /u(A) = for some 
quaternionic polynom,ial /i(A). Let B = PAP^^ be a similar matrix, with P 
a real matrix. Then fJ,{B) = 0. 

Proof. Let /x(A) = qiXq2^ ■ ■ ■ qk^qu+i be a monomial. Then /i(-B) = P fi{A)P^^ . 

D 

Notice that the same result is true when //(A) is a rational function. 

By permuting rows and columns (see Remark 13. 2p we deduce: 

Corollary 6.5. Let A 6e a 3 x 3 quaternionic matrix with some zero entry 
outside the diagonal. Then there exists a polynomial characteristic function 
/i such that fi{A) = 0. 



Example 6.6. Let us consider the matrix A 



"1 


i 


i 


i 


J 


k 





-1 


i_ 



It is real similar 



to 



IS 



j 


-1 





k 


j 


i 


i 


i 


1 



whose characteristic function is given by formula ([7]), that 



/x(A) = i + i{j - A) + (1 - X)i {k + (j - A)2) . 
Then the following equation holds: 

-AiA^ + AiAj + AkA + iA^ - iAj + A{i + j) - {i + k)A + {k - j)I = 0. 

6.3 Case n = 3,c ^ 

When c 7^ 0, the characteristic function of the matrix ^ is a rational function 
with a pole. We shall extend it to a map in the space of matrices in the 
following natural way. 

Let Xq = g — hc~^b be the pole of A. Let 



fo = f -he (a - Ao) 



qo = q- {r - Ao)c ^b. 
Lemma 6.7. The matrix XqI — A is invertible if and only if fo,qo 7^ 0. 

Proof By Corollary H^l Sdet(Ao/ - A) = \c\\qofo\. D 

Definition 6.8. We define jj: M{n,M.) — > M{n,M) as follows (see Defini- 
tion Ol): 



1. if XqI — B is invertible, then IJ,{B) = Qofol'^ 

2. otherwise, 

fi{B) = {Xol - B) [{pi - {r - B)c-\al - B)) - 

{ql _ (^/ _ B)c-H) (Ao - B)-^ [fl - hc-\al - B))] . 



The following Proposition completes the proof of Theorem A. 

Proposition 6.9. The map p in Def. \6.^ satisfies Hamilton- Cayley theo- 
rem, that is fJ,{A) = 



Proof. If XqI — ^ is not invertible, then ^{A) = qofol = by Lemma IGTTI 
Otherwise it suffices to prove that 



pi - {rl - A)c-\al - A) 



(8) 



equals 



{ql - {rl - A)c-^b) (Ao/ - A)-^ {fl - hc-^{al - A)) . (9) 

Lemma 6.10. A direct computation shows that the first term ^ is 

-ber\f -q+{r-a)c-^b + bc-^a-g) -bc'^h 

+ (r — g)c^^f + hc^^p p — fc^^b — hc^^q — (r — g)c^^{a — g) — / + [r — g)c^^h + hc^^{a — r) 
~qc^^ f —pc^^b + qc^^{a — g) —qc~^h 

We now want to compute the term ([9]). 

We start by computing {XqI — j4)~^ by Gaussian elimination. 

Let 

1 0" 

1 

c-i(Ao-a) 1 



10 



and 



Then 



where 



1 


0" 





1 





-c-^b 1 



(Ao/ - A)PiP2 



-c 

-/o -h 
-P* -qo Ao - r 



(10) 



p* =p- [Xq- r)c ^(Ao - a). 
The inverse of the matrix {\qI — A)PiP2 in (|10p can be computed by hand; 



it is 

where 

It follows that 



B 



fo'hc-^ 



— 1 * — 1 * J" — 1 — 1 

Qo ri Qo p /o -go 
-c-i 



n''=P*fo'hc-'-{Xo-r)c-\ 
(Ao/ - A)-^ = P1P2B = 



-1 * 








% 



Moreover 



F = fI-hc-^{aI-A) 

while 

Q = ql-{rl - A)c-^b -- 



f hc-^b h 

hc-^f f-hc-^{a-g) hc-^h 

hc~^p hc~^q f — hc~^{a — r) 



q — {r — a)c^^b bc^^b b 

fc~^b q — (r — g)c''^b hc^^b 

pc^^b qc~^b q 



We have to compute (HD, that is QP1P2BF. 

First we compute {PiP2B)F. For instance, its first column is given by 



[{PiP2B)F]^ 



-qo^n*f + q^Vfo'hc-'f - q^^hcr^p 
+c-^bq^'n*f - c-^f - c-^bq^'p*f^'hc-^f + cr^bq^^hC^p 



11 



Now we check for instance the entry (1, 1) of the matrix Q{PiP2BF). 
We have 



[QiPiP2BF)]\ = 

bc-'b {-q^'n*f + qoVfo'hc-'f - q^^hc-^p) + 

b {+c-hq^^n*f - c-'f - c-'bq.Vf^^hc-'f + c-'bq^'hc~'p) = 

-bc~'f 



which indeed is the entry (1, 1) in Corollary 16.101 

The other entries are computed in a similar way. D 

(1 ^ -J\ 
Example 6.11. Let A = I i —I k . The pole is Aq = —2 and 

V 1 -1 J / 

/i(Ao) = — 5 + 8j. For A ^ 2, the characteristic function is 
^(A) = -(2 + A) (2 + A(-l +i) - AjA + (-1 + i - AA:)(2 + A)-ii(2 - A)) . 
With the notations of the of proof of Proposition 16.91 it is 

/ -3 3i 

[XqI - A)-^ = (1/12) I 2i-j-k -8 + 2i + j + 3k 2 + 2i + 4k 
\l + i-j -3-i + 2j + k 2-2 + j + k 

/ -j l-i + 3k 1 
P=i -k 3-i-3j -i 
\ -k -2j + k i 

-1 + i — k j i 

Q= \ j -l + i + k I 

-k k -1 



and 



i 1 k 
F=\ I 3i j 

-i i 2i — k 



We have P - Q{\qI - A)-^F = 0. 



12 



7 Final remarks 

In this Section we discuss a different approach to the definition of charac- 
teristic functions for left eigenvalues. 

In order to clarify concepts, let us briefly comment the same problem 
but for right eigenvalues. Let c{A) G A4(2n,C) be the complex form of the 
matrix A G A^(n,IHI) (see Section [2|). Then, as it is well known, the right 
eigenvalues of A are the quaternions qzq~^, where q gM, q ^ 0, and 2; is a 
complex eigenvalue of c{A). It follows: 

Theorem 7.1 ([9]). Let p{z) = det{c{A) - zl) = Y.kLo^^kz'', Ck G M, be the 
characteristic polynomial of c{A). Then p{A) = X^^^q'^'s^'^ ~ *-*■ 

Now, let A = x+jy, with x, y G C, be a left eigenvalue oi A. Equivalently, 
the matrix c{A — XI) is not invertible. It follows that the left eigenvalues 
are the roots of the function a: C x C ^> M given by 



a{x,y) = det 



X -xl -Y + yl 
Y -yl X-xI 



ill) 



Let A = X + jY, with X,Y £ M{n,C)- Then Hamilton-Cayley theorem 
could be stated as o'{X, Y) = 0, provided this has a meaning. However we 
have the following counterexample even for n = 2. 



Example 7.2. Let A 



i 
J 



Let X = xi + ix2, y = yi + iy2- Then 



a{x,y) = l + {xi+X2 + yi+y2) -4x22/1. 



On the other hand, it is X 



i 




and Y 



1 




Yl 





1 



and Y2 = 0, then a{X,Y) 





1 
-3 

1 



, so Xl 



0, X2 
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